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Smoke  Plume  Trajectory  from  In  Situ  Burning  of  Crude  Oil  in  Alaska 


K.B.  McGrattan,  A.D.  Putorti,  W.H.  Twilley,  D.D.  Evans 
Building  and  Fire  Research  Laboratory 
National  Institute  of  Standards  and  Technology 
U.S.  Department  of  Commerce 
Gaithersburg,  MD  20899 


Abstract 

Experimentation,  analysis,  and  modeling  have  been  performed  to  predict  the 
downwind  dispersion  of  smoke  resulting  from  in  situ  burning  of  oil  spills.  North  Slope  and 
Cook  Inlet  crude  oils  are  burned  on  water  in  a 1.2  meter  diameter  pan.  Smoke  yields 
were  found  to  be  11 .6%  ±1 .0  for  North  Slope  crude,  and  9.2%  ±0.6  for  Cook  Inlet  crude 
by  mass  of  fuel  consumed,  with  a 95%  confidence  interval.  Burning  rates  and  smoke 
aerosol  size  distributions  are  also  measured,  and  found  similar  to  previous  work  with 
different  crude  oils.  Derivation  of  scaling  factors  for  predicting  the  burning  rates  and 
smoke  yields  of  large  scale  fires  are  guided  by  previous  experiments  with  Louisiana  crude 
oil.  Scaled  burning  rates  and  smoke  yields  are  supplied  as  input  parameters  for  the  LES 
(Large  Eddy  Simulation)  model,  version  2.0,  of  windblown  smoke  transport  over  flat 
terrain.  For  weather  conditions  appropriate  for  the  Cook  Inlet  and  North  Slope  areas, 
model  results  are  presented  which  predict  downwind  dispersion  and  ground  level 
concentrations  of  the  fire  generated  particulate  matter.  The  model  predicts  that  ground 
level  particulate  concentrations  in  excess  of  150  ^g/m^  are  limited  to  strips  5 km  long  and 
1 km  wide  downwind  of  the  fire  for  all  meteorological  conditions  considered. 

Key  words:  crude  oil;  oilspills;  pool  fires;  smoke;  fire  plumes 
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EXECUTIVE  SUMMARY 


A principal  concern  in  the  decision  to  use  in  situ  burning  as  an  oil  spill  mitigation 
technique  is  the  anticipated  trajectory  of  the  smoke  plume  and  the  dispersion  of  the  par- 
ticulate matter  generated  by  the  fire.  In  this  report  we  describe  the  Large  Eddy  Simulation 
(LES)  model  of  smoke  transport,  version  2.0,  which  may  be  used  to  predict  the  concen- 
tration of  combustion  products  downwind  of  the  fire.  Experiments  have  been  conducted 
to  determine  the  heat  release  rate  and  smoke  yield  from  two  types  of  Alaskan  crude 
oils.  North  Slope  and  Cook  Inlet,  burning  in  a 1.2  meter  diameter  pan.  Based  on  results 
from  mesoscale  experiments  performed  in  a 15  meter  square  pan  at  the  US  Coast  Guard 
mesoscale  burn  facility  in  Mobile  Bay,  Alabama,  the  laboratory  results  for  the  Alaskan 
crude  oils  are  extrapolated  to  mesoscale  to  be  used  as  input  for  the  LES  plume  trajectory 
model. 

The  mathematical  model  consists  of  solving  the  conservation  equations  of  mass, 
momentum  and  energy  which  describe  the  steady-state,  convective  transport  of  heated 
gases  introduced  into  the  atmosphere  by  a continuously  burning  fire  and  blown  by  a 
uniform  ambient  wind.  We  do  not  model  the  fire  itself,  but  rather  the  plume  of  smoke 
which  emanates  from  it.  Only  the  convective  heat  release  rate  and  smoke  yield  per  unit 
area  of  burning  surface  are  required  from  experiments.  Then,  as  input  to  the  model,  one 
need  only  specify  the  area  of  burning  surface,  the  wind  speed  and  its  fluctuations,  and 
the  temperature  stratification  of  the  atmosphere.  We  consider  two  burn  areas,  232  and 
465  m^  (2500  and  5000  ft^),  which  are  roughly  semi-circular  and  contained  in  the  apex  of 
a boom.  It  is  estimated  that  burns  covering  these  areas  will  consume  about  13,000  and 
27,000  gallons  of  oil  per  hour  and  generate  500  and  1 ,000  megawatts  of  heat,  respectively. 
For  the  Cook  Inlet  crude,  about  9%  of  the  mass  of  fuel  consumed  will  be  converted  into 
particulate  matter;  for  the  North  Slope  crude,  roughly  12%.  Wind  speeds  in  the  range 
from  about  5 to  25  knots  are  considered.  The  temperature  profiles  used  in  the  calculations 
correspond  to  typical  summer  and  winter  conditions  of  the  Cook  Inlet  and  North  Slope 
regions. 

In  processing  the  results  of  the  model,  special  attention  was  given  to  the  down- 
wind and  lateral  extent  of  ground  level  particulate  concentrations  in  excess  of  150 
averaged  over  one  hour.  Although  there  is  no  formal  guidance  available  about  safe  lev- 
els of  particulate  from  oil  fires,  150  fig/m^  averaged  over  24  hours  is  the  upper  level 
established  under  air  quality  control  regulations  in  Alaska  [13].  For  the  range  of  parame- 
ters discussed  above,  the  calculations  show  that  particulate  concentrations  found  at  the 
ground  downwind  of  oil  spill  burns  will  not  exceed  150  fig/m^  beyond  5 km,  nor  outside 
of  a strip  approximately  1 km  wide  along  the  centerline  of  the  plume  trajectory.  Given  the 
assumptions  and  uncertainties  inherent  in  the  model,  it  is  reasonable  to  assign  a factor  of 
safety  of  2 to  these  distances.  This  is  especially  prudent  since  the  model  does  not  at  the 
moment  consider  terrain  effects  or  coastal  weather  conditions  which  could  possibly  lead 
to  higher  ground  level  concentrations  than  those  predicted  by  the  model. 
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PART  1:  EXPERIMENT  AND  ANALYSIS 


1.1  Background 

This  study  is  part  of  a continuing  effort  at  NIST  aimed  at  understanding  and 
evaluating  in  situ  burning  as  a marine  oil  spill  response  tool.  Extensive  studies  of  crude 
oil  layers  burning  on  water  have  been  conducted  at  various  scales  ranging  from  0.085 
meters  to  1 7.2  meters  in  diameter,  with  several  different  crude  oils  [1 ,2].  During  previous 
studies,  however,  burns  greater  than  2.0  meters  in  diameter  have  been  conducted 
exclusively  with  Louisiana  crude  oil.  Since  full  scale  data  is  not  available  for  Cook  Inlet 
and  North  Slope  crude  oil  burns,  we  assume  that  the  behavior  of  these  crude  oils  is 
similar  to  that  of  Louisiana  crude  oil,  which  at  this  time  is  the  only  oil  studied  extensively 
at  a large  range  of  fire  sizes. 

In  order  to  model  the  burning  and  smoke  plume  formation  resulting  from  the  use 
of  in  situ  burning,  it  is  important  to  know  the  burning  rate  of  the  fuel.  The  burning  rates 
of  fires  of  different  scales  may  be  compared  if  the  burning  rate  of  fuel  is  expressed  as  the 
mass  loss  rate  per  unit  area  of  fuel  surface,  ie.  kg/(s»m^).  The  regression  rate  of  the  fuel 
surface  is  also  useful  for  comparison  purposes,  and  is  related  to  the  burning  rate  per  unit 
area  via  oil  density. 

The  fraction  of  fuel  mass  consumed  by  the  fire  that  is  emitted  as  smoke  aerosol, 
defined  as  smoke  yield,  is  another  parameter  necessary  for  the  modeling  of  smoke 
plumes.  In  this  investigation,  the  terms  smoke  aerosol,  smoke,  and  smoke  particulate  are 
all  synonymous  and  are  defined  as  the  condensed  phase  components  of  the  combustion 
products  [3]. 

Uncertainties  stated  in  this  report  represent  a 95%  confidence  interval  [4,5]. 

1.2  Experimental  Apparatus 

The  experimental  series  was  conducted  at  NIST’s  large  scale  fire  experimentation 
facility.  A steel  pan  measuring  1 .2  meters  in  diameter  and  0.30  meters  deep  was  placed 
under  the  large  calorimeter  hood.  The  bottom  of  the  pan  is  equipped  with  a pressure 
transducer  used  to  monitor  the  mass  of  oil  and  water  in  the  pan.  At  the  start  of  each 
burn,  a 25  ±4  mm  thick  layer  of  oil  was  added  on  top  of  245  ±4  mm  of  water,  placing  the 
top  surface  of  the  oil  approximately  30  mm  from  the  lip  of  the  pan.  The  experimental 
series  consisted  of  1 Louisiana,  5 North  Slope,  and  4 Cook  Inlet  crude  oil  burns. 

The  calorimeter  is  equipped  for  oxygen  depletion  calorimetry,  which  is  based  on 
the  principle  that  most  burning  materials  emit  13.1  MJ  of  heat  per  kg  of  oxygen 
consumed.  [6]  The  mass  of  oxygen  consumed  is  derived  from  the  decrease  in  oxygen 
concentration  in  the  combustion  products  flowing  through  the  exhaust  stack  as  compared 
to  ambient  air.  The  drop  in  oxygen  concentration  is  multiplied  by  the  volume  flow  rate  of 
exhaust  gases  through  the  stack,  and  the  density  of  oxygen  at  the  exhaust  temperature, 
resulting  in  the  rate  of  oxygen  consumption.  This  capability  allowed  the  heat  release  rate 
of  the  fire  to  be  recorded  at  regular  intervals  during  the  experiments  independently  of  fuel 
consumption  rate. 

An  elevated  platform  was  constructed  adjacent  to  the  calorimeter  to  physically 
support  smoke  sampling  equipment.  Sampling  probes  were  inserted  into  the  duct,  their 
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lengths  kept  short  to  avoid  heat  losses  and  smoke  aerosol  agglomeration.  Flow  rates 
were  selected  to  provide  for  isokinetic  sampling  within  the  exhaust  duct  in  order  to 
prevent  a bias  toward  small  or  large  smoke  particle  sizes.  Instrumentation  was  provided 
to  measure  smoke  particle  size  distribution  and  smoke  particle  mass.  Grab  samples  were 
collected  for  gas  chromatographic  analysis  of  CO  and  COg  levels  in  the  exhaust  for 
smoke  yield  calculations. 

1.3  Crude  Oil  Analysis 

Samples  of  each  of  the  crude  oils  supplied  by  the  Alaskan  Department  of 
Environmental  Conservation  (ADEC)  were  sent  to  a commercial  laboratory  for  elemental 
analysis.  The  laboratory  provided  the  mass  fractions  of  carbon,  hydrogen,  nitrogen,  and 
sulfur  in  the  oils.  The  results  are  in  Table  1 . 


Table  1.  Elemental  analysis  of  crude  oils. 


Element 

Cook  Inlet 
(mass  fraction) 

North  Slope 
(mass  fraction) 

Louisiana  [1] 
(mass  fraction) 

Carbon 

0.8432  ±0.0058 

0.8511  ±0.0058 

0.862  ±0.002 

Hydrogen 

0.1247  ±0.0045 

0.1210  ±0.0045 

0.134  ±0.016 

Nitrogen 

<0.0050 

<0.0050 

— 

Sulfur 

0.0019  ±0.0034 

0.0127  ±0.0034 

0.000  ±0.040 

The  density  of  each  of  the  oils  was  determined  at  NIST  by  the  mechanical 
oscillator  method.  This  technique  derives  the  density  of  a constant  volume  sample  by 
measuring  the  period  of  electromagnetically  induced  vibrations.  [7]  Densities  and  API 
gravities  of  the  oils,  along  with  viscosities,  and  pour  points  where  available  [8]  are  listed 
in  Table  2. 

1.4  Qualitative  Burning  Characteristics 

As  illustrated  in  previous  studies  with  other  oils  [9],  the  crude  oils  exhibited  various 
burning  stages  when  burned  in  the  1 .2  meter  pan.  These  stages  include:  ignition,  growth, 
steady  burning,  boiling,  and  extinction.  Cook  Inlet  and  North  Slope  crude  oils  also 
followed  this  sequence,  and  all  of  the  experiments  were  similar  during  the  steady  burning 
phase,  in  which  measurements  of  burning  rate  and  smoke  yield  were  recorded. 

Differences  in  ignition  behavior,  however,  were  noted  between  the  two  different 
oils.  Cook  Inlet  crude,  for  example,  was  ignited  in  1 to  2 seconds  with  a propane  torch, 
resulting  in  rapid  flame  spread  over  the  surface  of  the  oil.  The  North  Slope  oil,  in 
contrast,  was  more  difficult  to  ignite,  requiring  the  torch  to  be  held  on  the  surface  of  the 
oil  for  20  seconds  in  some  burns.  The  flame  spread  rate  was  also  slower  for  North  Slope 
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Table  2.  Crude  oil  properties. 


Cook  Inlet 

North  Slope 

Louisiana 

Density  (kg/m^) 
at  26‘’C 

841  ±1 

876  ±1 

842  ±1  [2] 

API  Gravity 
(15/15“C) 

35.0 

28.5 

36.0 

Kinematic  Viscosity 
(mm^/s)  at  38  °C 

— 

16.0 

4.3 

Pour  Point  ("C) 

— 

-8 

-9 

crude,  with  some  portions  of  the  surface  "flashing"  repeatedly  before  sustaining  flame. 
The  ignition  behavior  of  North  Slope  crude  oil  appeared  to  be  a strong  function  of 
temperature  in  the  ambient  temperature  range  of  the  burn  series  which  varied  from  20 
to  25  degrees  Celsius. 

1.5  Heat  Release  Rates  and  Burning  Rates 

The  heat  release  rates  of  the  experiments  are  plotted  as  a function  of  time  in 
Figures  1 and  2.  These  two  graphs  illustrate  the  progression  of  the  fires  from  ignition  to 
extinction.  Ignition  takes  place  at  time  zero,  at  which  time  the  fire  enters  the  growing 
phase  indicated  by  an  increasing  heat  release  rate  over  time.  The  plateau  on  the  graphs 
is  indicative  of  the  steady  burning  phase  of  the  fire,  during  which  all  of  the  sampling  takes 
place.  Immediately  following  the  steady  phase  is  the  boiling  phase  of  the  fire,  the  start 
of  which  is  represented  by  the  sudden  increase  and  peaking  of  the  heat  release  rate. 
The  process  of  extinction  begins  when  the  heat  release  rate  decreases  from  the  peak 
heat  release  rate. 

Using  the  pressure  transducer  at  the  base  of  the  pan,  the  weight  of  crude  oil  and 
water  may  be  derived  at  any  time  during  the  experiment.  Converted  to  layer  thicknesses 
via  initial  oil  density,  the  change  in  the  oil  layer  thickness  may  also  be  found  as  a function 
of  time.  Oil  surface  regression  rates  and  mass  loss  rates  during  the  steady  burning 
regime  may  be  found  in  Table  3,  along  with  heat  release  rates  from  oxygen  consumption 
calorimetry. 

As  was  the  case  with  previous  research,  the  mass  loss  rates  were  found  to  be 
relatively  constant  over  the  entire  burn,  although  a measurable  increase  does  take  place 
during  the  boiling  phase.  [2] 
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Heat  Release  Rate  (kW) 


Figure  1.  Heat  release  rate  for  North  Slope  crude  oil  burns. 
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Figure  2.  Heat  release  rate  for  Cook  Inlet  crude  oil  burns. 
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Table  3.  Average  burning  rate  and  heat  release  rate,  steady  burning 
phase. 


Burn  no. 

Oil 

kg/(s»m^) 

kW/m^ 

kW 

mm/s 

1 

Louisiana 

0.027 

1100 

1300 

0.032 

2 

North  Slope 

0.032 

— 

— 

0.036 

3 

North  Slope 

0.032 

1240 

1400 

0.036 

4 

North  Slope 

0.029 

1060 

1200 

0.033 

5 

North  Slope 

0.032 

1150 

1300 

0.036 

6 

North  Slope 

0.026 

1150 

1300 

0.030 

7 

Cook  Inlet 

0.031 

1330 

1500 

0.037 

8 

Cook  Inlet 

0.037 

1330 

1500 

0.044 

9 

Cook  Inlet 

0.031 

1240 

1400 

0.037 

10 

Cook  Inlet 

0.032 

1240 

1400 

0.038 

average* 

Louisiana 

0.033 

±0.006 

1400 

±200 

1600 

±400 

0.039 

±0.009 

average 

North  Slope 

0.030 

±0.005 

1150 

±70 

1300 

±80 

0.034 

±0.006 

average 

Cook  Inlet 

0.033 

±0.005 

1280 

±50 

1450 

±60 

0.039 

±0.006 

Values  are  the  average  of  burn  1 and  a previous  experiment  [2]  with  values 
of  0.039  kg/(s»m^),  1630  kW/m^,  1800  kW,  0.046  mm/s. 


1 .6  Smoke  Yield 

The  quantity  of  smoke  emitted  by  the  combustion  process  is  an  important 
parameter.  Smoke  emissions  may  be  expressed  as  smoke  yield,  which  is  defined  as  the 
mass  of  smoke  aerosol  generated  per  mass  of  fuel  consumed.  Various  methods  may  be 
used  to  measure  smoke  yield  including  the  flux  method,  the  carbon  balance  method,  and 
the  specific  light  extinction  method.  [10]  In  small  scale  experiments,  the  results  of  the 
three  methods  were  within  6%.  [1]  Since  the  carbon  balance  method  is  the  only  one  of 
the  three  methods  that  can  be  used  in  both  the  laboratory  and  in  the  field,  it  is  the 
method  used  in  this  study. 
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The  carbon  balance  method  is  based  on  the  assumption  that  the  mass  of  carbon 
in  the  fuel  consumed  is  equal  to  the  mass  of  carbon  in  the  airborne  combustion  products. 
Given  that  the  smoke  aerosols  are  predominately  carbon,  it  may  be  assumed  that  the 
carbon  consumed  by  the  fire  is  converted  to  particulate,  CO,  and  CO2  present  in  the 
smoke  plume.  [2]  Isokinetic  probes  in  the  exhaust  hood  stack  are  used  to  sample  a 
quantity  of  smoke,  which  is  passed  through  a filter  which  removes  the  smoke  aerosols. 
In  order  to  average  out  small  transients  in  burning  behavior,  the  sample  is  taken  over  a 
period  of  4 minutes  during  the  steady  phase  of  the  fire.  Due  to  the  experimental 
procedure  of  sampling  combustion  products  from  the  exhaust  stack,  it  is  necessary  that 
the  smoke  aerosols  and  combustion  gases  be  well  mixed  within  the  exhaust  stack,  and 
the  smoke  be  sampled  over  a time  great  enough  to  average  out  fire  fluctuations.  These 
additional  assumptions  are  valid  for  the  experimental  conditions,  and  the  smoke  yield  is 
found  from  the  following  expression,  as  derived  from  previous  work  [2]: 


smoke  yield  = 


m,  + 12  n (Axco,  + Axco) 


(1) 


where: 


n= 


AXco2- 


Axco- 


mass  of  smoke  collected  on  the  filter 

mass  fraction  of  carbon  in  the  fuel 

moles  of  gas  flowing  through  the  filter 

difference  in  volume  fraction  of  CO2  between  the 

sample  and  the  background 

difference  in  volume  fraction  of  CO  between  the 
sample  and  the  background 


Smoke  yields  for  Louisiana,  Cook  Inlet,  and  North  Slope  crude  oils  have  been 
calculated  for  each  experiment  and  listed  in  Table  4. 


1.7  Particle  Size  Distribution 

The  particle  size  distribution  of  smoke  aerosols  is  important  for  the  evaluation  of 
smoke  plume  dispersion  and  health  effects.  Due  to  the  irregular  shape  of  smoke 
particles,  which  consist  of  chains  of  small  spherules  [9],  the  particles  are  typically 
classified  by  any  one  of  a number  of  different  equivalent  diameters.  In  this  study,  the 
aerodynamic  properties  of  the  smoke  aerosol  are  of  primary  concern,  necessitating  the 
use  of  equivalent  aerodynamic  diameter.  The  equivalent  aerodynamic  diameter  of  an 
irregularly  shaped  particle  is  the  diameter  of  a smooth  spherical  particle,  having  a unit 
density  of  1 g/cm^,  with  the  same  terminal  velocity  as  the  smoke  particle  falling  in  air 
under  the  influence  of  gravity. 

Particle  size  was  measured  during  the  steady  portion  of  the  fires  using  two  particle 
impactors,  each  with  8 stages  and  a backup  filter.  The  impactors  are  built  in  such  a way 
that  a number  of  overlapping  perforated  disks,  or  stages,  force  the  particles  to  change 
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Table  4.  Average  smoke  yield,  steady  burning  phase. 


Burn  No. 

Oil 

Smoke  Yield 

1 

Louisiana 

0.101 

2 

North  Slope 

0.102 

4 

North  Slope 

0.119 

5 

North  Slope 

0.114 

6 

North  Slope 

0.115 

7 

Cook  Inlet 

0.088 

8 

Cook  Inlet 

0.093 

9 

Cook  Inlet 

0.093 

10 

Cook  Inlet 

0.095 

average* 

Louisiana 

0.105  ±0.010 

average 

North  Slope 

0.116  ±0.010 

average 

Cook  Inlet 

0.092  ±0.006 

*Average  of  burn  1 and  a previous  experiment  [2]  with  a value  of  0.109. 


direction  and  follow  the  bulk  flow  of  air  through  the  holes  in  the  stage.  Some  of  the 
particles,  however,  are  unable  to  navigate  this  path  due  to  their  momentum  and  are 
deposited  on  the  stage.  Each  stage  is  designed  to  collect  particles  of  a certain  size 
range,  specified  by  the  manufacturer  as  the  outpoint  diameter  of  the  stage,  which  is  the 
aerodynamic  diameter  that  is  collected  at  50%  efficiency.  The  cutpoint  diameter 
decreases  with  each  successive  stage  of  the  impactor.  The  cutpoint  sizes  for  the 
impactors  are  listed  in  Table  5. 

Two  different  types  of  impactor  were  used  in  this  experimental  series.  The  first 
was  a full  size  unit,  operating  at  a flow  rate  of  28.4  liters  per  minute.  Due  to  the  large 
size  of  this  impactor,  the  mass  of  each  stage  is  greater,  as  is  the  quantity  of  particulate 
collected.  This  combination  of  features  helps  decrease  the  uncertainty  in  the  particle  size 
measurements. 

The  second  impactor  used  in  this  study  was  primarily  designed  to  be  worn  by 
workers  to  monitor  particulate  exposure  levels  in  the  workplace.  Due  to  the  small  size 
of  this  impactor,  it  has  also  been  used  as  part  of  an  airborne  smoke  plume  sampling 
package  for  large  scale  fires.  As  a consequence  of  the  size  and  weight  advantages. 
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Table  5.  Impactor  stage  outpoint  sizes. 


Large  Impactor 
(28.4  liters/minute) 

Small  Impactor 
(2.0  liters/minute) 

Stage  1 (^im) 

9.0 

21.3 

Stage  2 (^im) 

5.8 

14.8 

Stage  3 (^m) 

4.7 

9.8 

Stage  4 (p,m) 

3.3 

6.0 

Stage  5 (^im) 

2.1 

3.5 

Stage  6 (^m) 

1.1 

1.55 

Stage  7 (^im) 

0.65 

0.93 

Stage  8 (^im) 

0.43 

0.52 

Back-up  Filter  (pm) 

0 

0 

however,  the  flow  rate  of  air  sampled  by  the  impactor  is  only  2.0  liters  per  minute.  The 
capacities  of  the  stages  are  also  small,  limiting  the  quantities  of  particulate  that  may  be 
collected. 

The  smoke  particle  sizes  for  the  oils  have  been  plotted  as  cumulative  distributions 
in  Figures  3 through  6.  The  distributions  are  very  repeatable,  especially  with  the  high  flow 
rate  impactor.  Notice  that  both  the  North  Slope  and  the  Cook  Inlet  oils  yield  almost 
identical  results  as  the  Louisiana  burn.  In  reference  to  health  concerns,  the  plots  indicate 
that  approximately  80%  to  99%  of  the  smoke  aerosols  collected  have  aerodynamic 
diameters  of  10  jim  (microns)  or  less.  This  is  significant  since  the  quantity  of  particulate 
matter  below  10  ^im  cut  size  is  a parameter  used  by  the  United  States  Environmental 
Protection  Agency  to  gauge  particulate  pollution  in  ambient  air.  [11] 
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Smoke  aerosol  size  distribution  for  North  Slope  crude  oil  burns  using 
the  high  flow  rate  impactor. 


12 


Cumulative  Mass  (%) 


99.9 


Figure  4. 
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Smoke  aerosol  size  distribution  for  Cook  Inlet  crude  oil  burns  using 
the  high  flow  rate  impactor. 
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Figure  5. 
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Smoke  aerosol  size  distribution  for  North  Slope  crude  oil  burns  using 
the  personal  impactor. 
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Figure  6. 
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Smoke  aerosol  size  distribution  for  Cook  Inlet  crude  oil  burns  using 
the  personal  impactor. 
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PART  2:  DETERMINATION  OF  SMOKE  PLUME  CALCULATION  INPUTS 


2.0  Fire  Scenario 

We  must  now  formulate  scenarios  of  interest  to  which  the  smoke  dispersion  model 
will  be  applied.  It  is  assumed  that  crude  oil  from  a spill  is  being  burned  within  the  roughly 
semi-circular  apex  area  of  a boom.  We  will  consider  two  burn  areas,  232  m^  (2500  fr) 
and  465  m^  (5000  ft^). 

2.1  Derivation  of  Scaling  Factors 

Comparing  the  burning  rates  of  Cook  Inlet  and  North  Slope  crude  oils  to  previous 
work  with  Louisiana  crude  oil  is  necessary  to  predict  the  large  scale  behavior  of  the 
former  oils.  Analyses  of  the  Louisiana  burn  conducted  in  this  study,  and  of  experiments 
conducted  previously  have  been  graphed  in  Figure  7.  This  plot  clearly  demonstrates  an 
increase  in  burning  rate  as  the  diameter  of  the  fire  increases.  At  fire  diameters  larger 
than  6.88  meters,  however,  the  burning  rate  plateaus,  and  becomes  independent  of  fire 
size.  From  this  graph,  a scaling  factor  may  be  defined  relating  the  laboratory  burning  rate 
to  the  large  scale  burning  rate,  ie. 

large  scale  burning  rate  = scaling  factor  x laboratory  burning  rate 

Due  to  the  similarities  of  the  three  crude  oils  at  burns  of  1 .2  meters  in  diameter,  the 
scaling  factor  is  assumed  to  be  the  same  for  all  of  the  oils.  The  burning  rate  of  the  three 
crude  oils  is  thereby  expected  to  increase  by  a factor  of  1 .7  when  the  fire  diameter  is 
increased  from  1 .2  meters  to  diameters  greater  than  6.88  meters,  which  would  include 
both  fires  stated  in  the  fire  scenario.  Since  heat  release  rates  from  field  experiments  are 
derived  from  the  burning  rate,  the  heat  release  rate  will  have  the  same  scaling  factor  as 
burning  rate.  The  results  of  applying  the  scaling  factor  of  1.7  to  burning  rate,  and 
consequently  to  heat  release  rate,  are  listed  in  Table  6. 

Table  6.  Estimated  large  scale  burning  rate,  heat  reiease  rate,  and  smoke 
yield. 


Oil 

Burning  rate 

Heat  release  rate 

Smoke  yield 

Cook  Inlet 

0.056  kg/(s»m^) 

5.9  gal/(hr*ft^) 

2180  kW/m^ 
6.90x10^  Btu/(hr«ft2) 

0.092 

North  Slope 

0.051  kg/(s»m^) 

5.1  gal/(hr*ft^) 

1960  kW/m^ 
6.20x10^  Btu/(hr*ft^) 

0.116 
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Figure  7.  Louisiana  crude  oil  surface  regression  rate.  Previous  [2]  and  new 
data. 

A scaling  factor  may  also  be  derived  for  smoke  yield.  Smoke  yield  has  been 
plotted  as  a function  of  fire  diameter  in  Figure  8.  From  this  graph,  it  can  be  seen  that  the 
smoke  yield  is  relatively  constant  for  fire  diameters  greater  than  1 .2  meters.  Therefore, 
the  scaling  factor  would  be  1 .0,  and  the  smoke  yields  for  the  large  scale  fire  scenarios 
would  be  expected  to  be  the  same  as  the  smoke  yields  from  the  laboratory  burns.  The 
scaling  factor  for  smoke  yield  has  been  applied  to  the  laboratory  data,  and  the  results 
listed  in  Table  6. 
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Figure  8.  Louisiana  crude  oil  smoke  yield.  Previous  [2]  and  new  data. 

Smoke  particle  size  measurements  from  all  three  of  the  oils  showed  nearly 
identical  results,  indicating  that  the  trends  observed  for  the  Louisiana  burns  may  be 
extrapolated  to  both  Cook  Inlet  and  North  Slope  oils.  The  particle  distributions  from  the 
1 .2  meter  burns  are  also  comparable  to  previous  results  for  17.2  meter  burns,  illustrated 
by  Figure  9,  indicating  a scaling  factor  of  1 .0.  Therefore,  a similar  distribution  would  be 
expected  for  a full  scale  burn  of  North  Slope  or  Cook  Inlet  crude  oils,  resulting  in  80%  to 
99%  of  the  smoke  aerosols  meeting  the  under  10  p.m  criterion.  This  distribution 
represents  the  smoke  aerosol  in  the  portion  of  the  plume  close  to  the  fire,  however,  and 
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Figure  9.  Smoke  aerosol  size  distribution  for  Louisiana  crude  oil  using  the 
personal  impactor.  [2] 

would  be  expected  to  change  as  the  plume  progressed  and  the  mechanisms  of 
agglomeration  and  settling  became  significant. 

2.2  Plume  Model  Inputs 

The  general  scaled  results  from  the  previous  section  are  now  used  to  derive  the 
specifics  of  the  fire  scenarios.  As  an  example,  the  first  case  would  involve  a Cook  Inlet 
crude  oil  burn  of  465  m^.  Based  on  experimental  data  we  assume  that  the  heat  and 
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mass  release  from  the  fire  is  more  or  less  steady.  Given  the  heat  release  and  burning 
rate  values  from  Table  6,  and  the  burning  area  of  the  fire,  the  total  heat  release  rate  and 
burning  rate  for  the  fire  are  calculated  in  Equations  3 and  4 respectively. 

Q = 465  X 2^80  kW/m^  = WOO  MW  (^) 

M = 465/772  0.056 /cg/(s-m2)  = 20.0 kg! s (4) 

The  "heat  loading"  of  the  fire,  or  the  convective  heat  release  rate,  is  defined  as  the  total 
heat  release  of  the  fire  minus  the  fraction  released  as  radiation.  Limited  radiation  data 
from  large  fires  [12]  suggest  a radiative  fraction  of  approximately  10%,  resulting  in  the 
following  calculation  of  heat  loading; 

heat  loading  = WOO  MW  x 0.90  = 900  MW  (5) 

The  smoke  yield  of  the  burn  is  expected  to  be  9.2%  from  Table  6,  so  the  smoke 
particulate  mass  flux  rate,  referred  to  as  "mass  loading",  is  given  by: 

mass  loading  = 0.092  x 20.0  kg! s = 2.39  kgis  (®) 

This  procedure  is  repeated  for  the  other  three  scenarios,  with  the  results  recorded  in 
Table  7. 
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Table  7.  Fire  Scenario  Summary 


Scenario  1 
Cook  Inlet 
Crude 

Scenario  2 
Cook  Inlet 
Crude 

Scenario  3 
North  Slope 
Crude 

Scenario  4 
North  Slope 
Crude 

Fire 

Area 

465 

(5000  ft^) 

232 

(2500  ft^) 

465 

(5000  ft^) 

232  m^ 
(2500  ft^) 

Fire 

Output 

Power 

1000  MW 
(3.45x10^ 
Btu/hr) 

500  MW 
(1.73x10® 
Btu/hr) 

910  MW 
(3.10x10® 
Btu/hr) 

450  MW 
(1.55x10® 
Btu/hr) 

Heat 

Loading 

900  MW 
(3.11x10® 
Btu/hr) 

450  MW 
(1.56x10® 
Btu/hr) 

820  MW 
(2.79x10® 
Btu/hr) 

405  MW 
(1.40x10® 
Btu/hr) 

Burning 

Rate 

26.0  kg/s 
(2.94x10"^ 
gal/hr) 

13.0  kg/s 
(1.47x10"^ 
gal/hr) 

23.7  kg/s 
(2.57x10"^ 
gal/hr) 

11.8  kg/s 
(1.29x10"^ 
gal/hr) 

Mass 

Loading 

2.39  kg/s 
(5.28  Ibm/s) 

1 .20  kg/s 
(2.63  Ibm/s) 

2.75  kg/s 
(6.07  Ibm/s) 

1 .37  kg/s 
(3.03  Ibm/s) 
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PART  3:  SMOKE  PLUME  TRAJECTORY  MODELING 


The  LES  (Large  Eddy  Simulation)  smoke  plume  trajectory  model  is  being  developed 
to  address  the  need  for  accurate  predictions  of  downwind  smoke  dispersion  from  large 
buoyant  sources  such  as  oil  spill  fires  to  support  the  decision  process  about  the  application 
of  in  situ  burning  in  response  to  oil  spills.  This  model  differs  from  most  of  the  atmospheric 
dispersion  models  in  use  today  because  it  is  a deferm/n/sf/c  rather  than  an  emp/r/ca/ model; 
that  is,  the  approach  taken  Is  to  solve  the  governing  equations  of  motion  directly  rather 
than  relying  on  empirical  formulae  which  approximate  the  extent  of  the  dispersion.  These 
empirical  models  typically  assume  the  pollutant  of  interest  to  be  Gaussian  distributed  in 
the  plane  perpendicular  to  the  direction  of  the  prevailing  wind.  The  parameters  defining 
the  distribution  are  estimated  from  experiments.  Unfortunately,  the  problem  of  in  situ 
burning  of  crude  oil  is  inappropriate  for  these  types  of  models  for  two  reasons:  (1)  The 
nature  of  the  “source”  is  different  than  what  is  normally  assumed,  a smokestack,  and  (2) 
the  size  of  the  source  is  well  beyond  those  considered  in  industrial  applications  and  thus 
outside  of  the  experimental  parameter  range. 

Because  the  LES  model  approaches  the  dispersion  problem  from  first  principles,  it 
is  more  flexible,  and  ultimately  easier  to  apply  because  it  reduces  the  number  of  empirical 
parameters  demanded  of  the  user.  The  long  range  goal  of  the  LES  modeling  study  is  to 
provide  high  resolution  predictions  that  can  account  for  the  actual  terrain  and  atmospheric 
conditions  in  the  vicinity  of  the  spill.  LES  version  2.0  used  in  this  study  is  limited  to 
predictions  of  downwind  smoke  concentration  for  uniform  winds  and  flat  terrain.  This 
version  represents  an  improvement  over  the  previous  version  [1]  with  the  inclusion  of 
atmospheric  dispersion,  more  efficient  computation,  and  clearer  presentation  of  results. 

3.1  Mathematical  Model 

We  assume  that  the  plume  may  be  described  in  terms  of  the  steady-state  convective 
transport  by  a uniform  ambient  wind  of  heated  gases  and  particulate  matter  introduced 
into  the  atmosphere  by  a continuously  burning  fire.  We  do  not  Intend  to  model  the  fire 
itself,  but  rather  the  plume  of  smoke  which  emanates  from  it.  Only  the  convective  heat 
release  rate  and  the  fraction  of  the  fuel  converted  to  particulate  matter  need  be  specified 
from  experiments.  The  simulation  begins  several  fire  diameters  downwind  of  the  fire, 
where  the  flow  field  is  characterized  by  relatively  small  temperature  perturbations,  minimal 
radiation  effects,  and  a velocity  field  dominated  by  the  prevailing  wind.  In  this  region  the 
plume  gases  ascend  to  a point  in  the  atmosphere  of  neutral  buoyancy,  and  then  gradually 
disperse.  The  trajectory  of  the  plume  is  governed  by  the  ambient  wind,  the  atmospheric 
stratification  and  the  level  of  turbulent  motion.  Since  it  is  not  our  objective  to  calculate  the 
local  meteorology,  it  is  assumed  that  this  information  is  available.  The  wind  is  taken  as 
uniform,  but  the  temperature  may  vary  with  height  according  to  a prescribed  profile. 

Given  these  assumptions,  the  mathematical  model  of  a smoke  plume  consists  of 
the  conservation  equations  of  mass,  momentum  and  energy  which  govern  the  absolute 
temperature  T,  pressure  p,  density  p,  and  crosswind  velocity  components  (v,  w)  in  a 
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plane  {y,z)  normal  to  the  direction  (x)  of  the  uniform  ambient  wind  (See  Figure  10).  It  is 


Figure  10:  Three  dimensional  view  of  a computed  smoke  plume  in  the  first  few 
kilometers  of  its  development.  The  height  of  the  viewbox  is  1 kilometer,  the  length 
8 kilometers,  and  the  crosswind  length  4 km.  The  wind  speed  is  6 m/s.  The  com- 
putation is  initialized  by  prescribing  the  temperature  and  particulate  distribution  in 
the  plane  spanned  by  the  y and  z coordinates.  Then  the  plume  is  constructed  as 
the  initial  plane  is  swept  downwind. 

convenient  to  divide  the  temperature  and  pressure  fields  into  mean  background  values 
Tq{z)  and  pq{z)  plus  perturbations  induced  by  the  fire,  f and  p,  respectively.  Similarly, 
the  density  p is  decomposed  into  an  ambient  density  po  and  a small  thermally  induced 
perturbation  p.  which  is  related  to  the  temperature  perturbation  through  the  equation  of 
state  taken  in  the  small  disturbance,  low  Mach  number  form  appropriate  to  this  problem 


T-To 

To 


P-  Po 
Po 


(7) 


The  ambient  density  is  related  to  the  background  pressure  through  the  hydrostatic  balance 


where  g represents  the  acceleration  of  gravity.  Assuming  that  the  perturbations  to  the 
background  temperature  and  density  are  small  beyond  a few  diameters  of  the  fire  down- 
wind of  the  firebed,  we  can  write  the  conservation  equations  describing  the  steady-state 
plume  in  the  Boussinesq  approximation  as  follows 

Conservation  of  mass 


dv  dw 


(9) 
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Conservation  of  lateral  and  vertical  momentum 


Po 

Po 


U-^  + v—  + w— 

^ ox  oy  oz 

\r^W  dw  dw 

U—  + v~  + w— 
ox  oy  oz 


d^v 


+ j;  + P9  = P 


Ur 


+ 


Conservation  of  energy 


far  dT  dT" 


dpQ  dTo  \ 

~ '’“‘^■57  j 


w = k 


d^f  d^f\ 

dy^  dz^  J 


(10) 

(11) 


(12) 


In  these  equations,  Cp  is  the  specific  heat  of  air,  k the  thermal  conductivity,  and  p the 
dynamic  viscosity.  The  uniform  ambient  wind  speed  U is  taken  to  be  constant  and  larger 
than  the  crosswind  velocity  components.  This  assumption  is  quite  realistic  several  flame 
lengths  downwind  of  the  firebed.  Since  U does  not  change,  there  is  no  need  for  a windward 
component  of  the  momentum  equations. 

The  three-dimensional,  steady  state  system  of  equations  above  can  be  considered 
as  a two-dimensional,  time-dependent  system  by  replacing  the  downwind  spatial  coordi- 
nate X with  a temporal  coordinate  t = xfLf.  We  then  have  an  initial  value  problem  in  which 
the  solution  is  initially  prescribed  in  a plane  perpendicular  to  the  direction  of  the  prevailing 
wind.  This  initial  plane  is  taken  to  be  a few  fire  diameters  downwind  of  the  fire.  Since  the 
details  of  the  firebed  are  not  being  simulated,  the  only  information  about  the  fire  required 
is  the  overall  convective  heat  release  rate  Q and  the  particulate  mass  flux  M.  The  initial 
temperature  and  particulate  distributions  in  the  plume  cross  section  are  assumed  to  be 
Guassian  and  satisfy  the  following  integrals 


roo  too  ^ 

/ / poCpUT  dzdy  — Q 

«/“00  Jq 

(13) 

too  too 

/ / ppU  dzdy  = M 

J —(XI  Jq 

(14) 

where  pp  is  the  concentration  of  particulate  matter,  defined  in  units  of  mass  per  unit 
volume.  In  the  simulations,  the  particulate  matter  is  represented  by  Lagrangian  particles 
which  are  advected  with  the  flow.  Initially,  the  crosswind  velocity  components  v and  w 
are  assumed  to  be  zero.  No-flux,  free-slip  boundary  conditions  are  prescribed  at  the 
ground,  consistent  with  the  assumed  uniformity  of  the  prevailing  wind  and  the  resolution 
limits  of  the  calculation.  At  the  outer  and  upper  edges  of  the  computational  domain,  the 
perturbation  temperature,  perturbation  pressure,  and  windward  component  of  vorticity  are 
set  to  zero.  Figure  10  shows  the  results  of  a sample  computation,  illustrating  the  position 
of  the  initial  slice  and  the  extent  of  the  computational  domain. 

Atmospheric  turbulence  effects  mixing  on  a wide  range  of  scales,  extending  to 
scales  which  are  smaller  than  the  resolution  of  the  calculations  performed  here.  There  are 
two  mechanisms  by  which  we  introduce  turbulent  atmospheric  motions  into  the  simulation. 
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Figure  11:  View  of  the  smoke  plume  from  underneath.  Notice  that  in  the  absence 
of  any  background  atmospheric  motion,  the  plume  is  characterized  by  two  large 
vortices  which  entrain  much  of  the  particulate  matter.  This  figure  shows  only  the 
fire-induced  flow  field. 

First,  the  small  scale  mixing  is  represented  by  a constant  eddy  viscosity  //,  which  is  taken 
to  be  three  orders  of  magnitude  greater  than  the  viscosity  of  air.  The  choice  is  governed 
by  the  desire  for  resolution  limits  in  the  five  to  fifteen  meter  range  which  are  needed  to 
capture  the  large  scale  fire  induced  eddy  motions.  This  requirement,  together  with  the 
knowledge  that  the  dissipative  effects  operate  at  a scale  Re~^  times  smaller  than  the 
overall  geometric  scale  (the  stabilization  height  of  the  plume  for  this  problem),  translates 
into  Reynolds  numbers  of  the  order  10^.  Thermal  conductivity  is  treated  in  a similar 
manner  to  viscosity,  thus  the  Prandtl  number  Pr  = Cpiijk  remains  of  order  unity.  The 
plume  depicted  in  Figure  10  represents  the  solution  of  Equations  (9)-(12)  with  these 
values  of  viscosity  and  thermal  conductivity. 

Unfortunately,  the  magnitude  of  the  viscosity  and  thermal  conductivity  included  in 
Equations  (9)-(12)  cannot  account  for  larger  scale  atmospheric  motions,  which  may  be 
thought  of  as  variations  of  the  prevailing  wind  over  a time  scale  of  minutes  to  hours. 
These  deviations  can  be  measured  with  an  anemometer,  and  are  introduced  into  the 
model  through  random  perturbations  to  the  trajectories  of  the  Lagrangian  particles  which 
represent  the  particulate  matter.  Thus,  the  motion  of  each  particle  is  governed  by  the  fire- 
induced  velocity  field  {U,v,w),  found  by  solving  the  conservation  equations  above,  plus 
a perturbation  velocity  field  (0,  v',  w')  which  represents  the  random  temporal  and  spatial 
variations  of  the  ambient  wind.  The  perturbation  velocity  components  are  derived  from 
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the  recursive  relations 


v\t  + 6t)  = R^i6t)v'{t)  + v"  (15) 

w'{t 6t)  = R.uj{6t)w\t)  + w"  (16) 

The  term  v"  is  a random  variable  with  Gaussian  distribution  whose  variance  is  that  of  v' 
multiplied  by  (1  -iilj),  ensuring  that  the  variance  of  the  velocity  does  not  change  from  one 
time  step  to  another.  The  term  w"  is  handled  the  same  way.  The  functions  Ry{6t)  and 
Ryj{6t)  are  Lagrangian  correlation  coefficients,  taken  as  exponentials 


R,{6t)  = e-**/"' 

(17) 

R„{St)  = e'*'/’’” 

(18) 

The  parameter  r is  indicative  of  the  period  of  atmospheric  fluctuations.  Appropriate  values 
for  various  meteorological  conditions  are  given  by  Draxler  [14].  Generally  speaking,  r is 
on  the  order  of  several  minutes,  and  it  does  not  have  much  of  an  effect  the  outcome  of 
the  simulations. 

The  initialization  of  the  recursive  relations  (15)  and  (16)  determines  the  overall 
level  of  turbulence  in  the  atmosphere.  The  initial  velocity  perturbations  i;'(0)  and  w'{0)  are 
randomly  assigned  from  a Gaussian  distribution  of  mean  zero  and  variance  aj  and  a^,  re- 
spectively. In  order  to  prescribe  these  initial  variations,  we  have  adopted  the  classification 
scheme  of  Pasquill  and  Gifford  [15]  to  categorize  local  meteorological  conditions  and  the 
level  of  atmospheric  motion.  These  classes  are  defined  in  Table  8.  Corresponding  to  each 

Table  8:  Key  to  Pasquill  stability  categories 


Surface 

Wind 

m/s 

Day 

Night 

Incoming  Solar  Radiation 

Cloud  Cover 

Strong 

Moderate 

Slight 

Low 

High 

<2 

A 

A-B 

B 

2-3 

A-B 

B 

C 

E 

F 

3-5 

B 

B-C 

C 

D 

E 

5-6 

C 

C-D 

D 

D 

D 

>6 

C 

D 

D 

D 

D 

Stability  Category  are  values  for  the  angular  variance  of  the  prevailing  wind  direction  in 
the  lateral  and  vertical  directions,  denoted  by  <t|  and  aj,  respectively.  The  initial  variance 
of  the  lateral  perturbation  velocity  t;'(0)  in  Equation  (15)  is  found  from  the  relation 

(Ty  = U sin  ae  (19) 

The  vertical  perturbation  velocity  w'  is  determined  in  a similar  way.  Reference  [1 6]  contains 
a discussion  of  these  parameters  and  methods  of  evaluating  them.  For  this  report,  we 
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have  adopted  the  following  values  of  oq  corresponding  to  the  Pasquill  Stability  Categories 
A-P;  25°,  20°,  15°,  10°,  5°,  and  2.5°,  respectively.  The  values  of  are  taken  to  be  16°, 
12°,  10°,  6°,  3°,  and  2°,  respectively.  These  values  of  wind  variation  vary  from  one  type 
of  terrain  to  another.  In  application,  best  results  can  be  expected  if  local  measurements 
of  average  wind  speed  and  variance  are  used  in  the  calculations.  These  values  generally 
vary  with  height.  Since  version  2.0  of  the  model  assumes  a constant  vertical  wind  profile, 
it  is  best  to  choose  the  values  which  correspond  to  the  expected  stabilization  height  of  the 
plume.  Also,  the  averaging  should  be  done  over  a time  period  which  roughly  corresponds 
to  the  expected  burning  time.  The  suggested  values  above  can  be  applied  to  burn  times 
of  ten  minutes  to  an  hour. 

3.2  Numerical  Method 

If  the  prevailing  wind  speed  is  greater  than  the  velocities  induced  in  the  crosswind 
direction  by  the  fire,  then  the  downwind  spatial  component  x may  be  considered  a tem- 
poral variable,  and  Equations  (9)-(12)  may  be  solved  with  a time  marching  procedure 
which  tracks  a planar  slice  of  particulate  matter  through  the  steady-state  plume.  The 
initial  temperature  distribution  of  the  plume  cross  section  is  assumed  to  be  Guassian;  its 
center  and  standard  deviation  are  estimated  from  observed  mesoscale  experiments  [1]. 
The  equations  are  solved  in  nondimensional  form  using  a relatively  simple  finite  difference 
technique.  All  spatial  derivatives  are  centrally  differenced,  and  the  solution  is  advanced  in 
time  (/.e.  the  initial  distributions  of  temperature  and  particulate  matter  are  tracked  down- 
wind) with  a second  order  Runge-Kutta  scheme.  The  computational  domain  spanning  the 
crosswind  plane  is  a rectangle  whose  ratio  of  vertical  to  lateral  lengths  is  usually  1 to  4. 
This  area  is  divided  into  rectangular  cells  (preferably  square)  with  a uniform  spacing  in 
each  direction.  For  the  calculations  reported  here,  a grid  with  64  cells  in  the  vertical  direc- 
tion and  256  cells  in  the  lateral  direction  was  used.  This  grid  resolves  the  crosswind  plane 
into  cells  from  5 to  20  meters  on  a side,  depending  on  the  scalings  applied.  On  current 
generation  work  stations,  this  is  a modest  calculation,  requiring  about  10  megabytes  of 
system  memory  and  10  to  20  minutes  of  CPU  time,  depending  on  the  specific  processor 
used. 
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PART  4:  RESULTS 


Part  2 of  this  report  contains  a discussion  of  the  derivation  of  input  parameters  for 
the  Large  Eddy  Simulation  model.  Briefly,  it  is  assumed  that  spilled  crude  oil  is  being 
burned  in  the  confines  of  a boom.  We  mainly  consider  two  fire  sizes,  one  corresponding 
to  a burning  area  of  465  m^  (5000  ft^)  and  the  other  232  m^  (2500  ft^).  (There  are  also 
a few  runs  with  a 10000  ft^  fire  for  comparison  purposes.)  We  assume  that  the  heat  and 
mass  release  rates  presented  in  Table  7 are  more  or  less  steady  for  burn  times  on  the 
order  of  one  hour. 

Figure  12  presents  the  results  from  a sample  run.  First,  Figure  12(A)  is  a typical 
temperature  profile  for  the  Cook  Inlet  region  in  the  winter.  The  effect  of  the  temperature 
inversion  on  the  plume  rise  is  evident  in  Figure  12(B)  which  shows  the  crosswind  extent 
of  hour-averaged  plume  particulate  concentrations  of  150  and  300  //g/m^  5 kilometers 
downwind  of  the  burn  site.  Figure  12(C)  is  a similar  plot  showing  the  downwind  extent  of 
these  concentration  levels  in  the  plume.  Figure  12(D)  displays  the  concentrations  at  the 
ground  level,  where  the  darker  shaded  areas  indicate  values  higher  than  150  /^g/m^.  The 
term  “ground  level”  refers  to  the  first  10  to  20  meters  of  the  atmosphere,  reflecting  the 
resolution  of  the  finite  difference  approximation  of  the  conservation  equations  described 
above.  Finally,  Figure  12(E)  is  included  to  quantify  the  previous  contour  plot.  It  shows  the 
decrease  in  the  ground  concentration  along  the  plume  centerline. 

For  each  region.  Cook  Inlet  or  North  Slope,  we  consider  the  burning  of  crude  oil 
from  that  particular  region,  i.e.  North  Slope  crude  oil  is  burned  on  the  North  Slope.  The 
heat  release  rates  and  the  smoke  yields  for  the  two  types  of  oil  considered  in  this  report  are 
similar  and  do  not  dramatically  change  the  results  of  the  computations.  For  example,  the 
case  presented  in  Figure  13  is  the  same  as  that  in  Figure  12  except  that  we  consider  the 
burning  of  North  Slope  crude  instead  of  Cook  Inlet  crude  in  the  vicinity  of  Cook  Inlet.  There 
is  very  little  change  in  the  final  result.  The  North  Slope  crude  has  a lower  heat  release 
rate  and  a higher  particulate  yield,  thus  its  burning  produces  a slightly  less  elevated  plume 
which  has  higher  particulate  concentrations  than  that  produced  from  the  burning  of  Cook 
Inlet  crude.  Given  the  assumptions  and  simplifications  inherent  in  the  model,  this  small 
difference  in  crude  oil  properties  is  of  little  consequence. 

Table  9 summarizes  the  results  of  the  simulations.  We  consider  four  different 
meteorological  conditions,  corresponding  to  the  Cook  Inlet  and  North  Slope  regions  in 
the  summer  and  winter.  Since  the  ultimate  objective  of  the  modeling  is  to  predict  the 
downwind  and  lateral  extent  of  harmful  concentrations  of  particulate  matter,  there  is  a 
column  entitled  “Extent  of  150  /^g/m^”  which  indicates  the  farthest  downwind  extent  of 
ground  level  concentrations  in  excess  of  150  ng/m^  averaged  over  one  hour.  Ground 
level  concentrations  in  excess  of  300  /xg/m^  did  not  extend  beyond  two  kilometers  for 
all  cases  shown;  and  most  often  not  even  beyond  one  kilometer.  These  concentrations 
represent  averaged  values  over  time  periods  on  the  order  one  hour.  The  Appendix 
includes  a synopsis  of  all  the  runs  listed  in  Table  9. 
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Finally,  note  that  Table  9 includes  a few  simulations  of  a larger  fire,  one  which 
would  cover  an  area  of  about  930  (10,000  ft^).  It  is  clear  that  even  though  the  fire 

size  is  double  the  area  of  the  465  case,  and  the  heat  and  mass  release  rates  are 
double,  the  particulate  concentrations  at  the  ground  are  not  double.  The  reason  for  this 
is  that  the  plume  will  rise  higher  when  the  fire  is  larger.  The  higher  elevation  of  the  plume 
causes  greater  dispersion,  and  thus  the  particulate  matter  is  spread  over  a larger  area. 
To  consider  fires  even  larger  than  the  ones  considered  here,  the  model  would  have  to  be 
altered  to  account  for  the  fact  that  multiple  plumes  might  arise  from  a fire  spread  over  a 
large  area. 
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Figure  12:  Sample  results  for  a typical  LES  plume  model  run.  This  case  considers 
the  burning  of  Cook  Inlet  crude  oil  in  the  vicinity  of  the  Cook  Inlet. 
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(B)  CROSSWINO  CONCENTRATION  AT  5000.  METERS 
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Figure  13:  Sample  results  for  a typical  LES  plume  model  run.  This  case  considers 
the  burning  of  North  Slope  crude  oil  in  the  vicinity  of  the  Cook  Inlet,  and  should  be 
compared  to  the  previous  example  in  which  Cook  Inlet  crude  is  burned  under  the 
same  conditions. 
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Table  9:  Summary  of  LES  smoke  plume  trajectory  model  runs 


Location 

Season 

Stability 

Class 

Wind 

Speed 

(m/s) 

Extent  of 
150 
(km) 

Burning  Area  of  232 

(2500  ft2) 

Cook  Inlet 

Summer 

C 

4 

<1 

Cook  Inlet 

Summer 

D 

8 

<1 

Cook  Inlet 

Summer 

D 

12 

1.5 

Cook  Inlet 

Winter 

C 

4 

<1 

Cook  Inlet 

Winter 

D 

8 

<1 

Cook  Inlet 

Winter 

D 

12 

2 

North  Slope 

Summer 

C 

4 

2.5 

North  Slope 

Summer 

D 

8 

3 

North  Slope 

Summer 

D 

12 

2.5 

North  Slope 

Winter 

C 

4 

3 

North  Slope 

Winter 

D 

8 

4 

North  Slope 

Winter 

D 

12 

2.5 

Burning  Area  of  465 

(5000  ft^) 

Cook  Inlet 

Summer 

C 

4 

<1 

Cook  Inlet 

Summer 

D 

8 

<1 

Cook  Inlet 

Summer 

D 

12 

1 

Cook  Inlet 

Winter 

C 

4 

<1 

Cook  Inlet 

Winter 

D 

8 

<1 

Cook  Inlet 

Winter 

D 

12 

2 

North  Slope 

Summer 

C 

4 

1 

North  Slope 

Summer 

D 

8 

4 

North  Slope 

Summer 

D 

12 

4 

North  Slope 

Winter 

C 

4 

5 

North  Slope 

Winter 

D 

8 

5 

North  Slope 

Winter 

D 

12 

4.5 

Burning  Area  of  930 

(lOOOOfL^) 

Cook  Inlet 

Summer 

D 

8 

<1 

Cook  Inlet 

Winter 

D 

8 

2 

North  Slope 

Summer 

D 

8 

2 

North  Slope 

Winter 

D 

8 

7 
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PART  5:  CONCLUSION 


North  Slope  and  Cook  Inlet  crude  oils  have  been  burned  on  water  in  a 1.2  meter 
pan  with  an  oil  layer  thickness  of  25  ±4  mm.  The  smoke  yields  of  these  fires  were  found 
to  be  11.6%  ±1.0  and  9.2%  ±0.6  per  mass  of  fuel  consumed, respectively.  The  burning 
rate  of  North  Slope  crude  was  0.030  ±0.005  kg/s-m^,  while  the  burning  rate  of  Cook  Inlet 
crude  oil  was  0.033  ±0.005  kg/s-m^.  Particles  with  equivalent  aerodynamic  diameters  of 

10  microns  or  less  were  found  to  represent  80  to  90%  of  the  total  mass  of  smoke  aerosol 
sampled.  It  should  be  noted  that  all  of  these  measurements  are  averages  over  the  steady 
burning  regime. 

Due  to  similarities  at  laboratory  scale,  the  scaling  factors  derived  by  comparing 
Louisiana  crude  oil  burns  of  various  sizes  are  applied  to  North  Slope  and  Cook  Inlet  crude 

011  burns.  The  scaling  factor  is  1 .7  for  burning  and  heat  release  rates,  and  1 .0  for  smoke 
yield  and  particle  size  distribution.  Therefore,  the  burning  rates  of  North  Slope  and  Cook 
Inlet  crude  oil  burns  at  equivalent  diameters  greater  than  6.88  meters  are  assumed  to 
be  0.051  kg/(s-m2)  and  0.056  kg/(s-m2),  respectively.  The  heat  release  rate  is  taken  as 
1960  kW/m^  for  North  Slope  and  2180  kW/m^  for  Cook  Inlet  crude  oils.  The  large  scale 
mass  based  smoke  yield  for  North  Slope  and  Cook  Inlet  crude  oils  are  assumed  to  be 
11.6%  and  9.2%,  respectively. 

We  use  these  extrapolated  values  of  heat  release  and  mass  flux  to  predict  the 
downwind  spread  of  combustion  products.  Despite  the  reasonably  high  resolution  of  the 
calculations,  uncertainties  need  to  be  recognized. 

• It  is  difficult  to  predict  the  convective  heat  release  rate  and  the  particulate  mass 
flux  from  the  fire  itself.  Even  though  we  have  confidence  in  the  laboratory  scale 
measurements,  there  is  some  uncertainty  in  extrapolating  the  results  to  the  scales 
anticipated  for  in  situ  burning  of  oil. 

• The  height  to  which  the  plume  rises  is  very  sensitive  to  the  background  temperature 
stratification.  The  model  results  vary  significantly  in  response  to  changes  in  the 
temperature  lapse  rate  of  a few  degrees  per  kilometer. 

• The  dispersion  parameters  reflecting  the  magnitude  of  the  wind  fluctuations  may 
vary  depending  on  the  terrain  of  the  burn  site.  The  values  used  for  the  simulations 
presented  here  are  based  on  empirical  estimates  of  plume  dispersion.  Better  results 
will  be  obtained  if  the  wind  fluctuations  are  measured  directly  with  the  use  of  an 
anemometer. 

• Finally,  there  is  very  little  experimental  data  with  which  to  verify  the  model. 

Nevertheless,  a study  of  a broad  range  of  conditions  has  yielded  useful  results.  After 
a number  of  simulations,  certain  patterns  and  trends  develop  which  instill  confidence  in 
the  overall  methodology.  For  the  conditions  considered  in  the  report  and  summarized  in 
Table  9,  we  found  that  the  LES  plume  trajectory  model  predicts  that  hour-averaged  ground 
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level  particulate  concentrations  of  150  ^g/m^  or  higher  do  not  extend  beyond  the  first  5 
kilometers  downwind  of  the  burn  site,  nor  do  these  levels  extend  outside  a path  of  about 
a kilometer  in  width.  Given  the  uncertainties  discussed  above,  these  estimates  can  be 
extended  by  a factor  of  safety  of  two  to  account  for  uncertainties  in  the  input  parameters. 

It  should  be  recognized  that  the  practical  application  of  in  situ  burning  of  oil  spills 
might  involve  more  than  one  fire  burning  at  a time,  with  varying  burn  times.  For  the 
moment,  the  LES  model  discussed  in  this  report  considers  only  the  burning  of  one  fire.  To 
account  for  multiple  burns  which  are  far  enough  away  from  each  other  to  be  considered 
independent,  it  is  reasonable  to  superimpose  the  results  of  several  simulations  to  yield  a 
total  concentration  at  any  location  as  the  sum  of  the  concentrations  from  the  respective 
runs.  If  the  fires  are  close  together  (in  relation  to  the  overall  size  of  the  plumes),  then  they 
may  be  considered  as  one  burn. 

Although  not  part  of  this  study,  another  practical  consideration  of  in  situ  burning  is 
the  time  required  to  clear  the  air  following  the  extinguishing  of  a burn.  This  question  is 
difficult  to  answer  because  of  terrain  effects.  It  is  reasonable  to  expect  the  air  above  30 
meters  in  altitude  to  clear  relatively  quickly,  and  the  speed  of  the  prevailing  wind  offers  a 
rough  guide  as  to  how  long  the  last  of  the  particulates  from  the  fire  pass.  However,  at  the 
surface,  where  knowledge  of  particulate  concentration  is  most  important,  it  is  difficult  to 
predict  with  the  present  version  of  the  model  how  long  it  will  take  for  the  air  to  clear.  The 
reason  for  the  difficulty  is  the  fact  that  if  the  terrain  is  rough,  or  if  it  is  an  urban  setting, 
pockets  of  high  concentration  may  get  “trapped”  between  buildings  or  in  a valley.  Local 
meteorological  and  topographical  conditions  would  thus  play  a crucial  role  in  determining 
the  time  required  for  the  air  to  clear. 
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APPENDIX:  LES  PLUME  TRAJECTORY  MODEL  RESULTS 


On  the  following  pages  we  present  the  results  of  simulations  of  possible  burns  off 
both  the  northern  and  southern  coasts  of  Alaska.  These  results  have  been  summarized  in 
Table  9.  The  temperature  profiles  were  provided  by  the  National  Weather  Service  office  in 
Alaska,  and  represent  typical  conditions  of  the  winter  and  summer  months.  For  each  case 
we  plot  the  vertical  temperature  profile,  along  with  plots  of  the  crosswind  and  downwind 
spread  of  the  particulate  matter.  We  focus  our  attention  on  particulate  concentrations  of 
150  and  300  //g/m^  averaged  over  one  hour.  A concentration  of  150  or  higher 
averaged  over  24  hours  is  considered  hazardous  to  human  health.  The  change  in  scale 
from  one  case  to  another  reflects  the  fact  that  the  wind  speed,  fire  size  and  temperature 
profile  determine  the  scale  heights  for  each  individual  case.  Thus,  it  is  important  to  note 
these  changes  in  scale  when  comparing  one  case  with  another.  A contour  plot  of  the 
ground  level  concentration  is  provided  to  indicate  the  downwind  and  lateral  spread  of  the 
particulate  matter  near  the  ground.  A plot  of  the  centerline  ground  concentration  is  also 
provided  to  aid  in  quantifying  the  contour  plot.  This  is  not  a plot  of  ground  deposition,  but 
rather  an  hour-averaged  airborne  concentration  of  particulate  matter  in  the  first  1 0 or  20 
meters  of  the  atmosphere. 

The  cases  shown  generally  correspond  to  weather  conditions  with  moderate  and 
slight  levels  of  daytime  incoming  solar  radiation.  (See  Table  8)  Generally  speaking,  for  the 
size  fires  that  we  are  considering,  wind  speeds  which  are  less  than  about  2 m/s  (4  knots) 
cannot  satisfy  the  assumption  that  the  prevailing  wind  be  stronger  the  induced  crosswinds. 
In  any  case,  we  have  found  that  a 6-1 0 m/s  (1 2-20  knots)  wind  in  combination  with  stability 
classes  C and  D will  produce  the  greatest  downwind  spread  of  particulate  matter,  and  thus 
we  focus  our  attention  on  this  regime.  Stability  classes  A and  B are  highly  dispersive  and 
generally  do  not  yield  high  ground  level  concentrations  beyond  about  a kilometer  of  the 
burn  site.  Classes  E and  F are  considered  night  time  conditions  and  are  inappropriate 
for  that  reason  and  because  the  wind  speeds  considered  are  too  high  for  their  use  (See 
Table  8). 


37 


DISTANCE  FROM  CENTERUNE  (METERS) 


COOK  INLET.  SUMMER 


WIND  SPEED  4.0  m/s 
STABIUTY  CLASS  C 


HEAT  LOADING  450.  MW 
MASS  LOADING  1.2  kg/s 


(A)  TEMPERATURE  PRORLE  (B)  CROSSWIND  CONCENTRATION  AT  5000.  METERS 


(C)  VERTICAL  CENTERUNE  CONCENTRATION 


(D)  GROUND  LEVEL  CONCENTRATION 


2500 

2000 

1500 

1000 

500 

0 

-500 

-1000 

-1500 

-2000 

-2500 


275.  250.  225.  200.  175.  150.  125.  100.  75.  50.  25.  0. 


2000 


4000 


6000 


BOOO  10000  12000 

DOWNWIND  DISTANCE  (METERS) 


14000 


16000 


18000 


20000 


2000 


4000 


6000 


ALL  CONCENTRATIONS  EXPRESSED  IN  UNITS  OF  /<g/m' 


8000  10000  12000 
DOWNWIND  DISTANCE  (METERS) 


14000 


16000 


18000 


20000 


Figure  14 


38 


DISTANCE  FROM  CENTERUNE  (METERS) 


COOK  INLET,  SUMMER 


WIND  SPEED  8.0  m/s 
STABIUTY  CLASS  D 


HEAT  LOADING  450.  MW 
MASS  LOADING  1.2  kg/s 


(A)  TEMPERATURE  PRORLE  (B)  CROSSWIND  CONCENTRATION  AT  5000.  METERS 


(C)  VERTICAL  CENTERUNE  CONCENTRATION 


2000 

1500 

1000 

500 

0 

-500 

-1000 

-1500 

-2000 


(0)  GROUND  LEVEL  CONCENTRATION 


0 2000  4000  6000  8000  10000  12000  14000  16000  18000  20000 

DOWNWIND  DISTANCE  (METERS) 


Figure  15 


39 


DISTANCE  FROM  CCNTERUNE  (METERS) 


COOK  INLET.  SUMMER 


WIND  SPEED  12.0  m/s 
STABIUTY  CLASS  D 


HEAT  LOADING  450.  MW 
MASS  LOADING  1.2  kg/s 


(C)  VERTICAL  CENTERUNE  CONCENTRATION 


1500 

1000 

500 

0 

-500 

-1000 

-1500 


(0)  GROUND  LEVEL  CONCENTRATION 


275.  250.  225.  200.  175.  150.  125.  100.  75.  50.  25.  0. 


0 2000  4000  6000  SOOO  10000  12000  14000  16000  1B000  20000 

DOWNWIND  DISTANCE  (METERS) 


0 2000  4000  6000  BOOO  10000  12000  14000  16000  16000  20000 

DOWNWIND  DISTANCE  (METERS) 

ALL  CONCENTRATIONS  EXPRESSED  IN  UNITS  OF 


Figure  16 


40 


DISTANCE  EROM  CENTERUNE  (METERS) 


COOK  INLET,  WINTER 


WIND  SPEED  4.0  m/s 
STABIUTY  CLASS  C 


HEAT  LOADING  450.  MW 
MASS  LOADING  1.2  kg/s 


(A)  TEMPERATURE  PROHLE 


TEMPERATURE  (C) 


(B)  CROSSWIND  CONCENTRATION  AT  SOOO.  METERS 


(C)  VERTICAL  CENTERUNE  CONCENTRATION 


1600 
1200 
800 
AOO 
0 

-400 
-800 
-1200 
-1600 

0 2000  4000  6000  8000  10000  12000  14000  16000  18000  20000 

DOWNWIND  DISTANCE  (METERS) 


(D)  GROUND  LEVEL  CONCENTRATION 


Figure  17 


41 


DISTANCE  FROM  CENTERUNE  (METERS) 


COOK  INLET,  WINTER 


WIND  SPEED  8.0  m/s 
STABIUTY  CLASS  D 


HEAT  LOADING  450.  MW 
MASS  LOADING  1.2  kg/s 


(A)  TEMPERATURE  PROHLE 


TEMPERATURE  (C) 


(B)  CROSSWIND  CONCENTRATION  AT  5000.  METERS 


(C)  VERTICAL  CENTERUNE  CONCENTRATION 


1200 
BOO 
400 
0 

-400 
-800 
-1200 

0 2000  4000  6000  8000  10000  12000  14000  16000  18000  20000 

DOWNWIND  DISTANCE  (METERS) 


(D)  GROUND  LEVEL  CONCENTRATION 


275.  250.  225.  200,  175.  150.  125.  100.  75.  50.  25.  0. 


0 2000  4000  6000  8000  10000  12000  14000  16000  18000  20000 

DOWNWIND  DISTANCE  (METERS) 

ALL  CONCENTRATIONS  EXPRESSED  IN  UNITS  OF  fiQ/rr? 


Figure  18 


42 


/ 


DISTANCE  FROM  CENTERUNE  (METERS) 


COOK  INLET,  WINTER 


WIND  SPEED  12.0  m/s 
STABIUTY  CLASS  D 


HEAT  LOADING  450.  MW 
MASS  LOADING  1.2  kg/s 


(A)  TEMPERATURE  PROHLE 


TEMPERATURE  (C) 


^500 

i/i 

u 400 

C 

5,300 

I 200 

UJ 

^ 100 
0 


(B)  CROSSWIND  CONCENTRATION  AT  5000.  METERS 


-1000  0 1000 
DISTANCE  FROM  CENTERUNE  (METERS) 


(C)  VERTICAL  CENTERUNE  CONCENTRATION 


1000 

800 

600 

400 

200 

0 

-200 

-400 

-600 

-800 

-1000 


(0)  GROUND  LEVEL  CONCENTRATION 


s’'" 

n 

III 

ill 

n 

275.  250.  225.  200.  175.  150.  125.  100.  75.  50.  25.  0. 


2000 


4000 


6000 


8000  10000  12000 
DOWNWIND  DISTANCE  (METERS) 


14000 


16000 


18000 


20000 


0 2000  4000  6000  8000  10000  12000  14000  16000  18000  20000 

DOWNWIND  DISTANCE  (METERS) 

AU  CONCENTRATIONS  EXPRESSED  IN  UNITS  OF 


Figure  19 


43 


DISTANCE  FROM  CENTERUNE  (METERS) 


NORTH  SLOPE.  SUMMER 


WIND  SPEED  4.0  m/s 
STABIUTY  CLASS  C 


HEAT  LOADING  410.  MW 
MASS  LOADING  1.4  kg/s 


(A)  TEMPERATURE  PROHLE 


(B)  CROSSWINO  CONCENTRATION  AT  5000.  METERS 


TEMPERATURE  (C) 


0 2000  4000  6000  BOOO  10000  12000  14000  16000  1B000  20000 

DOWNWIND  DISTANCE  (METERS) 


(D)  GROUND  LEVEL  CONCENTRATION 


1000 

800 

600 

400 

200 

0 

-200 

-400 

-600 

-BOO 

-1000 


275.  250.  225-  2Q0.W175.  75.  50.  25. 


2000 


4000 


6000 


BOOO  10000  12000 

DOWNWIND  DISTANCE  (METERS) 


14000 


16000 


1B000 


20000 


Figure  20 


44 


/ 


DtSTANCE  FROM  CENTERUNE  (METERS) 


NORTH  SLOPE,  SUMMER 


WIND  SPEED  8.0  m/s 
STABIUTY  CLASS  D 


HEAT  LOADING  410,  MW 
MASS  LOADING  1,4  kg/s 


(A)  TEMPERATURE  PROHLE 


(B)  CR0SSW1N0  CONCENTRATION  AT  5000.  METERS 


TEMPERATURE  (C) 


(C)  VERTICAL  CENTERUNE  CONCENTRATION 


(D)  GROUND  LEVEL  CONCENTRATION 


0 2000  4000  6000  BOOO  10000  12000  14000  16000  1BOOO  20000 

DOWNWIND  DISTANCE  (METERS) 


Figure  21 


45 


DISTANCE  FROM  CEMTERUNE  (METERS) 


NORTH  SLOPE.  SUMMER 


WIND  SPEED  12.0  m/s 
STABIUTY  CUSS  D 


HEAT  LOADING  410.  MW 
MASS  LOADING  1.4  kg/s 


(A)  TEMPERATURE  PROnUE 


TEMPERATURE  (C) 


(B)  CROSSWIND  CONCENTRATION  AT  SOOO.  METERS 


0 2000  4000  6000  BOOO  10000  12000  14000  16000  1B000  20000 

DOWNWIND  DISTANCE  (METERS) 


SOO 
600 
400 

200 
0 

-200 
-400 
-600 
-BOO 

0 2000  4000  6000  SOOO  10000  12000  14000  16000  16000  20000 

DOWNWIND  DISTANCE  (METERS) 


(0)  GROUND  LEVEL  CONCENTRATION 


■ ■ I 1 L. 


0 2000  4000  6000  6000  10000  12000  14000  16000  18000  20000 

DOWNWIND  DISTANCE  (METERS) 

ALL  CONCENTRATIONS  EXPRESSED  IN  UNITS  OF 


Figure  22 


46 


/ 


DISTANCE  FROM  CCNTERUNE  (METERS) 


NORTH  SLOPE,  WINTER 


WIND  SPEED  4.0  m/s 
STABIUTY  CLASS  C 


HEAT  LOADING  410.  MW 
MASS  LOADING  1.4  kg/s 


(A)  TEMPERATURE  PRORLE 


TEMPERATURE  (C) 


(B)  CROSSWIND  CONCENTRATION  AT  5000.  METERS 


(C)  VERTICAL  CENTERUNE  CONCENTRATION 


(D)  GROUND  LEVEL  CONCENTRATION 


0 2000  4000  6000  BOOO  10000  12000  14000  16000  18000  20000 

DOWNWIND  DISTANCE  (METERS) 


0 2000  4000  6000  8000  10000  12000  14000  16000  18000  20000 

DOWNWIND  DISTANCE  (METERS) 

ALL  CONCENTRATIONS  EXPRESSED  IN  UNITS  OF  yji/rT? 


Figure  23 


47 


DISTANCE  FROM  CENTERUNE  (METERS) 


NORTH  SLOPE,  WINTER 


WIND  SPEED  8.0  m/s 
STABILITY  CLASS  D 


HEAT  LOADING  410.  MW 
MASS  LOADING  1.4  kg/s 


(A)  TEMPERATURE  PRORLE 


(B)  CROSSWIND  CONCENTRATION  AT  5000.  METERS 


TEMPERATURE  (C) 


(C)  VERTICAL  CENTERUNE  CONCENTRATION 


1200 

1000 

800 

600 

400 

200 

0 

-200 

-400 

-600 

-BOO 

-1000 

-1200 


(0)  GROUND  LEVEL  CONCENTRATION 


0 2000  4000  6000  8000  10000  12000  14000  16000  1B000  20000 

DOWNWIND  DISTANCE  (METERS) 


0 2000  4000  6000  8000  10000  12000  14000  16000  18000  20000 

DOWNWIND  DISTANCE  (METERS) 

Aa  CONCENTRATIONS  EXPRESSED  IN  UNITS  OF  t^/rr? 


Figure  24 


48 


DISTANCE  FROM  CENTERUNE  (METERS) 


NORTH  SLOPE,  WINTER 


WIND  SPEED  12.0  m/s 
STABIUTY  CLASS  D 


HEAT  LOADING  410.  MW 
MASS  LOADING  1.4  kg/s 


(A)  TEMPERATURE  PROFILE 


TEMPERATURE  (C) 


(B)  CROSSWIND  CONCENTRATION  AT  5000.  METERS 


(C)  VERTICAL  CENTERUNE  CONCENTRATION 


(0)  GROUND  LEVEL  CONCENTRATION 


0 2000  4000  6000  8000  10000  12000  14000  16000  18000  20000 

DOWNWIND  DISTANCE  (METERS) 


0 2000  4000  6000  8000  10000  12000  14000  16000  18000  20000 

DOWNWIND  DISTANCE  (METERS) 

ALL  CONCENTRATIONS  EXPRESSED  IN  UNITS  OF  Atg/m’ 


Figure  25 


49 


DISTANCE  FROM  CENTERUNE  (MQERS) 


COOK  INLET,  SUMMER 


WIND  SPEED  4.0  m/s 
STABIUTY  CLASS  C 


HEAT  LOADING  900.  MW 
MASS  LOADING  2.4  kg/s 


(A)  TEMPERATURE  PRORLE 


(B)  CROSSWIND  CONCENTRATION  AT  5000.  METERS 


DISTANCE  FROM  CENTERUNE  (METERS) 


(C)  VERTICAL  CENTERUNE  CONCENTRATION 


3000 
2500 
2000 
1500 
1000 
500 
0 

-500 
-1000 
-1500 
-2000 
-2500 
-3000 

0 2000  4000  6000  8000  10000  12000  14000  16000  1B000  20000 

DOWNWIND  DISTANCE  (METERS) 


(0)  GROUND  LEVEL  CONCENTRATION 


0 2000  4000  6000  BOOO  10000  12000  14000  16000  18000  20000 

DOWNWIND  DISTANCE  (METERS) 

Aa  CONCENTRATIONS  EXPRESSED  IN  UNITS  OF  Mg/m’ 


Figure  26 


50 


DISTANCE  FROM  CENTERUNE  (METERS) 


COOK  INLET.  SUMMER 


WIND  SPEED  8.0  m/s 
STABIUTY  CLASS  D 


HEAT  LOADING  900.  MW 
MASS  LOADING  2.4  kg/s 


(A)  TEMPERATURE  PROnLX  (B)  CROSSWIND  CONCENTRATION  AT  SOOO.  METERS 


(C)  VERTICAL  CENTERUNE  CONCENTRATION 


2S00 

2000 

1500 

1000 

500 

0 

-500 

-1000 

-1500 

-2000 

-2500 


(0)  GROUND  LEVEL  CONCENTRATION 


275.  250.  225.  200.  175.  150.  125.  100.  75.  50.  25.  0. 


0 2000  4000  6000  8000  10000  12000  14000  16000  18000  20000 

DOWNWIND  DISTANCE  (METERS) 


0 2000  4000  6000  8000  10000  12000  14000  16000  18000  20000 

DOWNWIND  DISTANCE  (METERS) 

ALL  CONCENTRATIONS  EXPRESSED  IN  UNITS  OF  fig/rn^ 

Figure  27 


51 


DISTANCE  FROM  CENTERUNE  (METERS) 


COOK  INLET,  SUMMER 


WIND  SPEED  12.0  m/s 
STABIUTY  CLASS  D 


HEAT  LOADING  900.  MW 
MASS  LOADING  2.4  kg/s 


(A)  TEMPERATURE  PROFILE 


(B)  CROSSWIND  CONCENTRATION  AT  5000.  METERS 


(C)  VERTICAL  CENTERUNE  CONCENTRATION 


2000 

1500 

1000 

500 

0 

-500 

-1000 

-1500 

-2000 


(D)  GROUND  LEVEL  CONCENTRATION 


275.  250.  225.  200.  175.  150.  125.  100.  75.  50.  25.  0. 


0 2000  4000  6000  8000  10000  12000  14000  16000  18000  20000 

DOWNWIND  DISTANCE  (METERS) 


0 2000  4000  6000 

ALL  CONCENTRATIONS  EXPRESSED  IN  UNITS  OF 


8000  10000  12000 
DOWNWIND  DISTANCE  (METERS) 

Atg/m* 

Figure  28 


14000  16000 


18000  20000 


52 


/ 


DISTANCE  FROM  CENTERUNE  (METERS) 


COOK  INLET,  WINTER 


WIND  SPEED  4.0  m/s 
STABIUTY  CLASS  C 


HEAT  LOADING  900.  MW 
MASS  LOADING  2.4  kg/s 


(A)  TEMPERATURE  PROniX 


TEMPERATURE  (C) 


1000 

^ 800 
Ul 

y 600 

X 400 
o 

UJ 

z 200 
0 


(B)  CROSSWIND  CONCENTRATION  AT  SOOO.  METERS 


-2000  -1000  0 1000  2000 
DISTANCE  FROM  CENTERUNE  (METERS) 


(C)  VERTICAL  CENTERUNE  CONCENTRATION 


(D)  GROUND  LEVEL  CONCENTRATION 


2000 

1500 

1000 

500 

0 

-500 

-1000 

-1500 

-2000 


0 2000  4000  6000  8000  10000  12000  14000  16000  18000  20000 

DOWNWIND  DISTANCE  (METERS) 


0 2000  4000  6000  SOOO  10000  12000  14000  16000  18000  20000 

DOWNWIND  DISTANCE  (METERS) 

ALL  CONCENTRATIONS  EXPRESSED  IN  UNITS  OF  Mg/m’ 


Figure  29 


53 


DISTANCE  FROM  CCNTERUNE  (METERS) 


COOK  INLET,  WINTER 


WIND  SPEED  8.0  m/s 
STABIUTY  CLASS  D 


HEAT  LOADING  900.  MW 
MASS  LOADING  2.4  kg/s 


I 


(C)  VERTICAL  CENTERUNE  CONCENTRATION 


1600 

1200 

BOO 

400 

0 

-400 

-BOO 

-1200 

-1600 


(0)  GROUND  LEVEL  CONCENTRATION 


0 2000  4000  6000  8000  10000  12000  14000  16000  16000  20000 

DOWNWND  DISTANCE  (METERS) 


0 2000  4000  6000  BOOO  10000  12000  14000  16000  16000  20000 

DOWNWIND  DISTANCE  (METERS) 

Aa  CONCENTRATIONS  EXPRESSED  IN  UNITS  OF 


Figure  30 


54 


/ 


DISTANCE  FROM  CENTERUNE  (METERS) 


COOK  INLET,  WINTER 


WIND  SPEED  12.0  m/s 
STABIUTY  CLASS  D 


HEAT  LOADING  900.  MW 
MASS  LOADING  2.4  kg/s 


(A)  TEMPERATURE  PRORLE 


TEMPERATURE  (C) 


(B)  CROSSWIND  CONCENTRATION  AT  5000.  METERS 


(C)  VERTICAL  CENTERUNE  CONCENTRATION 


1200 

800 

400 

0 

-400 

-800 

-1200 


(D)  GROUND  LEVEL  CONCENTRATION 


0 2000  4000  6000  8000  10000  12000  14000  16000  18000  20000 

DOWNWIND  DISTANCE  (METERS) 


0 2000  4000  6000  8000  10000  12000  14000  16000  18000  20000 

DOWNWIND  DISTANCE  (METERS) 

Aa  CONCENTRATIONS  EXPRESSED  IN  UNITS  OF  MO/m’ 


Figure  31 


55 


DISTANCE  FROM  CENTERUNE  (METERS) 


NORTH  SLOPE,  SUMMER 


WIND  SPEED  4.0  m/s 
STABIUTY  CLASS  C 


HEAT  LOADING  820.  MW 
MASS  LOADING  2.8  kg/s 


(A)  TEMPERATURE  PROFILE 


TEMPERATURE  (C) 


(B)  CROSSWIND  CONCENTRATION  AT  5000.  METERS 


0 2000  4000  6000  BOOO  10000  12000  14000  16000  18000  20000 

DOWNWIND  DISTANCE  (METERS) 


(D)  GROUND  LEVEL  CONCENTRATION 


DOWNWIND  DISTANCE  (METERS) 


0 2000  4000  6000  8000  10000  12000  14000  16000  18000  20000 

DOWNWIND  DISTANCE  (METERS) 

Aa  CONCENTRATIONS  EXPRESSED  IN  UNITS  OF  Atg/m’ 

Figure  32 


56 


/ 


OfSTANCE  FROM  CCNTERUNE  (METERS) 


NORTH  SLOPE,  SUMMER 


WIND  SPEED  8.0  m/s 
STABIUTY  CLASS  D 


HEAT  LOADING  820.  MW 
MASS  LOADING  2.8  kg/s 


(A)  TEMPERATURE  PROFILE 


TEMPERATURE  (C) 


(B)  CROSSWIND  CONCENTRATION  AT  5000.  METERS 


0 2000  4000  6000  8000  10000  12000  14000  16000  1S000  20000 

DOWNWIND  DISTANCE  (METERS) 


1000 

800 

600 

400 

200 

0 

-200 

-400 

-600 

-800 

-1000 


(D)  GROUND  LEVEL  CONCENTRATION 


0 2000  4000  6000  8000  10000  12000  14000  16000  18000  20000 

DOWNWIND  DISTANCE  (METERS) 


0 2000  4000  6000  8000  10000  12000  14000  16000  18000  20000 

DOWNWIND  DISTANCE  (METERS) 

ALL  CONCENTRATIONS  EXPRESSED  IN  UNITS  OF  tii/rr? 


Figure  33 


57 


DISTANCE  FROM  CENTERUNE  (METERS) 


NORTH  SLOPE.  SUMMER 


WIND  SPEED  12.0  m/s 
STABIUTY  CLASS  D 


HEAT  LOADING  820.  MW 
MASS  LOADING  2.8  kg/s 


(A)  TEMPERATURE  PRORLE 


TEMPERATURE  (C) 


(B)  CROSSWIND  CONCENTRATION  AT  5000.  METERS 


DISTANCE  FROM  CENTERUNE  (METERS) 


(C)  VERTICAL  CENTERUNE  CONCENTRATION 


1000 

BOO 

600 

400 

200 

0 

-200 

-400 

-600 

-800 

-1000 


(D)  GROUND  LEVEL  CONCENTRATION 


0 2000  4000  6000  8000  10000  12000  14000  16000  18000  20000 

DOWNWIND  DISTANCE  (METERS) 


0 2000  4000  6000  8000  10000  12000  14000  16000  18000  20000 

DOWNWIND  DISTANCE  (METERS) 

ALL  CONCENTRATIONS  EXPRESSED  IN  UNITS  OF  MO/m’ 


Figure  34 


58 


/ 


DISTANCE  EROM  CENTERUNE  (METERS) 


NORTH  SLOPE,  WINTER 


WIND  SPEED  4.0  m/s 
STABIUTY  CLASS  C 


HEAT  LOADING  820.  MW 
MASS  LOADING  2.8  kg/s 


TEMPERATURE  (C) 


-2000  -1000  0 1000  2000 
DISTANCE  FROM  CENTERUNE  (METERS) 


2000 

1500 

1000 

500 

0 

-500 

-1000 

-1500 

-2000 


(D)  GROUND  LEVEL  CONCENTRATION 


0 2000  4000  6000  8000  10000  12000  14000  16000  1B000  20000 

DOWNWIND  DISTANCE  (METERS) 


Figure  35 


59 


DISTANCE  FROM  CENTERUNE  (METERS) 


NORTH  SLOPE.  WINTER 


WIND  SPEED  8.0  m/s 
STABIUTY  CLASS  D 


HEAT  LOADING  820.  MW 
MASS  LOADING  2.8  kg/s 


(A)  TEMPERATURE  PROHLE 


TEMPERATURE  (C) 


(B)  CROSSWIND  CONCENTRATION  AT  5000.  METERS 


(C)  VERTICAL  CENTERUNE  CONCENTRATION 


(0)  GROUND  LEVEL  CONCENTRATION 


0 2000  4000  6000  8000  10000  12000  14000  16000  1BOOO  20000 

DOWNWIND  DISTANCE  (METERS) 


Figure  36 


60 


DISTANCE  FROM  CENTERUNC  (METERS) 


NORTH  SLOPE.  WINTER 


WIND  SPEED  12.0  m/s 
STABIUTY  CLASS  D 


HEAT  LOADING  820.  MW 
MASS  LOADING  2.8  kg/s 


0 2000  4000  6000  BOOO  10000  12000  14000  16000  18000  20000 

DOWNWIND  DISTANCE  (METERS) 


1200 

800 

400 

0 

-400 

-BOO 

'1200 


(D)  GROUND  LEVEL  CONCENTRATION 


0 2000  4000  6000  6000  10000  12000  14000  16000  16000  20000 

DOWNWIND  DISTANCE  (METERS) 


0 2000  4000  6000  BOOO  10000  12000  14000  16000  18000  20000 

DOWNWIND  DISTANCE  (METERS) 

Aa  CONCENTRATIONS  EXPRESSED  IN  UNITS  OF 


Figure  37 


61 


DISTANCE  FROM  CENTERUNE  (METERS) 


COOK  INLET,  SUMMER 


WIND  SPEED  8.0  m/s 
STABIUTY  CLASS  D 


HEAT  LOADING  1800.  MW 
MASS  LOADING  4.8  kg/s 


(A)  TEMPERATURE  PROHLE 


(B)  CROSSWIND  CONCENTRATION  AT  5000.  METERS 


TEMPERATURE  (C) 


(C)  VERTICAL  CENTERUNE  CONCENTRATION 


2000 

2000 

1000 

0 

-1000 

-2000 

-3000 


(D)  GROUND  LEVEL  CONCENTRATION 


275.  250.  225.  200.  175.  150.  125.  100.  75.  50.  25.  0. 


0 2000  4000  6000  8000  10000  12000  14000  16000  18000  20000 

DOWNWIND  DISTANCE  (METERS) 


Figure  38 


62 


/ 


DISTANCE  FROM  CENTERUNE  (METERS) 


COOK  INLET,  WINTER 


WIND  SPEED  8.0  m/s 
STABIUTY  CLASS  D 


HEAT  LOADING  1800.  MW 
MASS  LOADING  4.8  kg/s 


(A)  TEMPERATURE  PRORLE 


TEMPERATURE  (C) 


(B)  CROSSWIND  CONCENTRATION  AT  5000.  METERS 


(C)  VERTICAL  CENTERUNE  CONCENTRATION 


(D)  GROUND  LEVEL  CONCENTRATION 


DOWNWIND  DISTANCE  (METERS) 


Figure  39 


63 


DISTANCE  FROM  CENTERUNE  (METERS) 


NORTH  SLOPE,  SUMMER 


WIND  SPEED  8.0  m/s 
STABIUTY  CLASS  D 


HEAT  LOADING  1640.  MW 
MASS  LOADING  5.6  kg/s 


(A)  TEMPERATURE  PRORLE 


TEMPERATURE  (C) 


-1000  0 1000 
DISTANCE  FROM  CENTERUNE  (METERS) 


(C)  VERTICAL  CENTERUNE  CONCENTRATION 


(D)  GROUND  LEVEL  CONCENTRATION 


0 2000  4000  6000  BOOO  10000  12000  14000  16000  1B000  20000 

DOWNWIND  DISTANCE  (METERS) 


Figure  40 


64 


DISTANCE  FROM  CENTERUNE  (METERS) 


NORTH  SLOPE,  WINTER 


WIND  SPEED  8.0  m/s 
STABIUTY  CLASS  D 


HEAT  LOADING  1640.  MW 
MASS  LOADING  5.6  kg/s 


TEMPERATURE  (C) 


-2000  -1000  0 1000  2000 
DISTANCE  FROM  CENTERUNE  (METERS) 


2000 

1500 

1000 

500 

0 

-500 

-1000 

-1500 

-2000 


(D)  GROUND  LEVEL  CONCENTRATION 


0 2000  4000  6000  8000  10000  12000  14000  16000  18000  20000 

DOWNWIND  DISTANCE  (METERS) 


Figure  41 


65 


/ 


N'ST-^  U.S.  DEPARTMENT  OF  COMMERCE 

(REV.  6-93)  NATIONAL  INSTITUTE  OF  STANDARDS  AND  TECHNOLOGY 

AOMAN  4.09 

MANUSCRIPT  REVIEW  AND  APPROVAL 

(ERB  USE  ONLY) 

ERB  CONTROL  NUMBER 

DIVISION 

PUBUCATION  REPORT  NUMBER 

NISTIR  5273 

CATEGORY  CODE 

INSTRUCTIONS:  ATTACH  ORIGINAL  OF  THIS  FORM  TO  ONE  (1)  COPY  OF  MANUSCRIPT  AND  SEND  TO 
THE  SECRETARY,  APPROPRIATE  EDITORIAL  REVIEW  BOARD. 

p06UCATION  date  NUMBER  PRINTED  PAGES 

October  1993 

TITLE  AND  SUBTITLE  (CITE  IN  FULL) 

Smoke  Plume  Trajectory  from  In  Situ  Burning  of  Crude  Oil  in  Alaska 


CONTRACT  OR  GRANT  NUMBER 


TYPE  OF  REPORT  AND/OR  PERIOD  COVERED 


AUTHOR(S)  (LAST  NAME,  FIRST  INITIAL,  SECOND  INITIAL) 

McGrattan,  Kevin  B.,  Putorti,  Anthony  D., 
Twilley,  William  H. , Evans,  David  D. 


PERFORMING  ORGANIZATION  (CHECK  (X)  ONE  BOX) 
NIST/GAITHERSBURG 
NIST/BOULDER 
JILA/BOULDER 


JL 


LABORATORY  AND  DIVISION  NAMES  (FIRST  NIST  AUTHOR  ONLY) 

Building  and  Fire  Research  Laboratory,  Fire  Safety  Engineering  Division 


SPONSORING  ORGANIZATION  NAME  AND  COMPLETE  ADDRESS  (STREET,  aTY,  STATE,  ZIP) 

Alaska  Dept,  of  Environmental  Conservation  (ADEC) 

410  Willowby  Avenue,  Suite  302 
Juneau,  AK  99801-1795 


PROPOSED  FOR  NIST  PUBUCATION 


JOURNAL  OF  RESEARCH  (NIST  JRES) 

MONOGRAPH  (NIST  MN) 

LETTER  CIRCULAR 

J.  PHYS.  & CHEM.  REF.  DATA  (JPCRD) 

NATL  STD.  REF.  DATA  SERIES  (NIST  NSRDS) 

BUILDING  SCIENCE  SERIES 

HANDBOOK  (NIST  HB) 

FEDERAL  INF.  PROCESS.  STDS.  (NIST  FIPS) 

PRODUCT  STANDARDS 

SPECIAL  PUBUCATION  (NIST  SP) 

UST  OF  PUBUCATIONS  (NIST  LP) 

OTHER 

TECHNICAL  NOTE  (NIST  TN) 

X 

NIST  INTERAGENCY/INTERNAL  REPORT  (NISTIR) 

PROPOSED  FOR  NON-NIST  PUBUCATION  (CITE  FULLY) 


U.S. 


FOREIGN 


PUaUSHING  MEDIUM 
PAPER 

DISKETTE  (SPECIFY) 
OTHER  (SPECIFY) 


□ 


CD-ROM 


SUPPLEMENTARY  NOTES 


ABSTRACT  (A  2000-CHARACTER  OR  LESS  FACTUAL  SUMMARY  OF  MOST  SIGNIFICANT  INFORMATION.  IF  DOCUMENT  INCLUDES  A SIGNIFICANT  BIBUOGRAPHY 
OR  UTERATURE  SURVEY,  CITE  IT  HERE.  SPELL  OUT  ACRONYMS  ON  FIRST  REFERENCE.)  (CONTINUE  ON  SEPARATE  PAGE.  IF  NECESSARV.I 

- Experimentation,  analysis,  and  modeling  have  been  performed  to  predict  the  downwind  dispersion  of  smoke  ' 
resulting  from  in  situ  burning  of  oil  spills.  North  Slope  and  Cook  Inlet  crude  oils  are  burned  on  water  in  a 1.2  meter 
diameter  pan.  Smoke  yields  were  found  to  be  11.6%  ±1.0  for  North  Slope  crude,  and  9.2%  ±0.6  for  Cook  Inlet  crude  by 
mass  of  fuel  consumed,  with  a 95%  confidence  interval.  Burning  rates  and  smoke  aerosol  size  distributions  are  also 
measured,  and  found  similar  to  previous  work  with  different  crude  oils.  Derivation  of  scaling  factors  for  predicting  the 
burning  rates  and  smoke  yields  of  large  scale  fires  are  guided  by  previous  experiments  with  Louisiana  crude  oil.  Scaled 
burning  rates  and  smoke  yields  are  supplied  as  input  parameters  for  the  LES  (Large  Eddy  Simulation)  model,  version  2.0, 
of  windblown  smoke  transport  over  flat  terrain.  For  weather  conditions  appropriate  for  the  Cook  Inlet  and  North  Slope  areas, 
model  results  are  presented  which  predict  downwind  dispersion  and  ground  level  concentrations  of  the  fire  generated  ' 
particulate  matter.  The  model  predicts  that  ground  level  particulate  concentrations  in  excess  of  150  pg/m^  are  limited  to 
strips  5 km  long  and  1 km  wide  downwind  of  the  fire  for  all  meteorological  conditions  considered. 


KEY  WORDS  (MAXIMUM  OF  9;  28  CHARACTERS  AND  SPACES  EACH;  SEPARATE  WITH  SEMICOLONS;  ALPH; 

Crude  oil;  fire  plumes;  oilspills;  pool  fires;  smoke 

VBETIC  ORDER;  CAPITAUZE  ONLY  PROPER  NAMES) 

AVAILABIUTY 

NOTE  TO  AUTHOR(S);  IF  YOU  DO  NOT  WISH  THIS 

X 

UNUMITED  1 1 FOR  OFFICIAL  DISTRIBUTION  - DO  NOT  RELEASE  TO  NTIS 

MANUSCRIPT  ANNOUNCED  BEFORE  PUBUCATION, 

ORDER  FROM  SUPERINTENDENT  OF  DOCUMENTS,  U.S.  GPO,  WASHINGTON,  DC  20402 

PLEASE  CHECK  HERt  | | 

_x_ 

ORDER  FROM  NTIS,  SPRINGFIELD,  VA  22161 

ELECTRONIC  FORM 


"1 

-i 

■I 

;1 


/ 


